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Abstract 


In this paper, a new modeling technique for fuel ejectors with high entrainment ratio, low pressure increment and over heated working gases in 
an anodic recirculation solid oxide fuel cell (SOFC) system is presented. By utilizing the thermodynamic, fluid dynamic principles and chemical 
constraints inside ejectors and employing a two-dimensional function to compute fluid velocity, the developed model involves no more than nine 
algebraic equations and this is very simple compared to all existing models. The detailed procedures for fuel ejector design and simulation are 
provided and its effectiveness is verified through simulation and compared with testing results. It shows that the proposed model is more accurate 
than presently available models, and therefore can be better used for ejector design and performance simulations. The ejector performances for 
both situations of stand-alone and integrated into the SOFC system are also studied. 


© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


Solid oxide fuel cell (SOFC) provides an effective means to 
convert chemical energy directly to electrical energy with high 
efficiency, flexible fuel utilization and very low pollutions [1—4], 
and this has drawn intensive research interests in recent years. 
A typical SOFC system includes a fuel cell stack and peripheral 
components such as fuel processor, heat exchanger, mixer and 
reformer. In the reformer, it is critical to have enough steam 
to prevent carbon deposition and provide sufficient heat for 
endothermic reforming reactions [5,6]. Since anodic exhausts 
are rich in steam and high in temperature, they can be recycled 
for fuel reforming by fans, blowers or ejectors. With no moving 
parts and less maintenance by using high pressure fuel gas as the 
primary fluid to suck the anodic exhausts, anodic recirculation 
using ejectors increases the SOFC system reliability compared 
with other schemes [7,8]. Since the cost of energy for fuel com- 
pression accounts for as high as 7% of the total cost of electricity 
[9], an extreme care should be taken in ejector design and opera- 
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tion for optimal system performance. Therefore, an accurate fuel 
ejector model for ejector design and the evaluation of on-design 
and off-design performances is essential. 

The behaviors of an ejector related to entrainment capac- 
ity, steam to carbon ratio (STCR) and outlet temperature are 
strongly influenced by the geometries and operational condi- 
tions of the ejector. The performances of ejectors have been 
studied for decades and several modeling methods are avail- 
able for ejectors design and performance evaluations [10-12]. 
However, most of the existing models are developed based 
on one-dimensional (1D) techniques for cooling and refriger- 
ation applications. These modeling methods will cause large 
errors to model fuel ejectors in SOFC systems due to the differ- 
ences in geometries, working fluids properties and operating 
conditions. Compared with traditional ejectors, fuel ejectors 
in an SOFC system are different mainly in the following 
aspects: 


e The diameter ratio of mixing chamber to nozzle throat is much 
bigger [11], due to the requirement of larger entrainment ratio. 

e The primary and secondary flows are overheated gases instead 
of saturated vapors. 
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Nomenclature 

A area (m2) 

Cp specific heat of gas at constant pressure 
(Jkg7! KT!) 

D diameter (m) 


FC fuel cell 


k specific heat ratio of gas 

m mass flow rate (kg s7!) 

M Mach number 

Mo molecular weight (kg mol~!) 
n molar flow rate (mol s7!) 

ny exponent in Eq. (12) 

P pressure (Pa) 

r,R radius (m) 

Rg gas constant (Jkg~! K7!) 

Ra universal gas constant (J mol`! K7!) 
T temperature (K) 


v, V velocity (m s7!) 


Greek symbols 

Wp isentropic coefficient of primary flow 

p density (kg m~?) 

w entrainment ratio, ms/mp 

exp coefficient accounting for friction loss during the 


mixing process 


Subscripts 

P primary flow (i.e. high pressure fuel) 

S secondary flow (i.e. low pressure anodic recycle 
gas) 

t nozzle throat 

0 ejector inlet 

1 primary flow at nozzle throat 

2 nozzle exit 

3 mixing chamber inlet 

4 diffuser exit 

Superscripts 

i chemical component 


e The secondary flow temperature is very high around 900°C. 
e The pressure increment of the secondary flow is much smaller. 


To the best of the authors’ knowledge, there are very few 
ejector models available for anodic recirculation SOFC systems. 
Marsano et al. developed an SOFC ejector model using 1D mod- 
eling technique for the evaluation of on-design and off-design 
performances [13]. Later, Ferrari et al. [14] improved the tech- 
nique by dividing the ejector into serials of calculation cells 
where the governing equations are numerically solved. How- 
ever, this model was also based on the 1D method. These 1D 
models assume that the velocity, pressure and temperature are 
uniformly distributed in the radial direction of the ejector. Based 
on this assumption 1D models can predict the performances of 


ejectors in traditional applications such as refrigeration within 
+15% [11]. However, since the flow area of the secondary flow 
in an anodic recirculation ejector is much greater than that of 
traditional application ejectors, the conventional 1D modeling 
techniques will cause larger errors when applied in modeling 
SOFC ejectors. 

In this paper, a simple yet accurate fuel ejector model for 
geometries design and performance evaluation is developed. The 
governing equations are derived through: (1) determining the 
primary mass flow rate by isentropic flow relations; (2) approx- 
imating the velocity distribution of secondary flow by a 2D 
function; (3) deriving a simple formula for secondary mass flow 
rate which is capable of dealing with the viscosity flow near the 
ejector inner wall and (4) establishing the energy conservation 
equation for the primary and secondary flows. The procedures 
for calculating ejector outputs are straightforward with no more 
than nine algebraic equations to be solved for both ejector design 
and performance simulation. To validate the effectiveness of 
the modeling technique for design, and performance evalua- 
tions at both on-design and off-design conditions, experimental 
numerical test and simulation studies have been carried out. Fur- 
thermore, the effects of pressure and temperature of primary flow 
and secondary flow, and chemical composition to the ejector per- 
formance on flow rate of anodic recycle gas, STCR, pressure, 
temperature and chemical composition at ejector outlet have also 
been studied. Finally, the performance of the fuel ejector in the 
SOFC system is investigated by varying the fuel inlet pressure, 
fuel cell pressure and fuel utilization. This model is expected 
to have wide applications in geometry design and performance 
evaluation of ejectors in anodic recirculation SOFC systems. 


2. Anodic recirculation SOFC system with fuel ejector 


A simplified sketch of a tubular SOFC is shown in Fig. 1 . It 
mainly consists of three components: an ejector, a reformer and 
a fuel cell stack. 

The working principle of the tubular SOFC can be briefly 
descried as follows: 


1. High pressure fuel (primary flow) passes through an ejector, 
where the low pressure anodic exhaust (secondary flow) is 
entrained and mixes with the primary flow. The resulting 


Exhausts 
pz 


Cathode 


Anode 


Fig. 1. Simplified sketch of a SOFC module. 
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mixing stream will shock in the diffuser to a higher pressure 

and then enter into the connected reformer. According to ref. 

[13], the function of the ejector is to re-circulate the anodic 

gas to 

e raise the secondary flow pressure to meet the FC pressure 
at the required level; 

e supply sufficient heat required for the reforming reactions 
in the reformer; 

e provide sufficient secondary flow rate to maintain a proper 
STCR to avoid carbon deposition in the reformer and FC 
stacks. 

2. Inside the reformer, highly endothermic reactions take place: 


Reforming: CH4+H20 < CO + 3H2 (la) 


Shifting: CO +H20 < CO? + H2 (1b) 


In the anodic recirculation SOFC system, the steam used 
for converting methane (or natural gas) into hydrogen in the 
reformer is supplied by the ejector. The required energy for 
reactions is provided by the sensible heat of the entrained 
high temperature secondary flow. 

3. The reformed fuel is fed to the anode side of FC stack, while 
air is supplied to the cathode side. In the cathode, oxygen 
ions passing through the electrolyte layer react with hydrogen 
diffusing through the anode to form steam, and the electrons 
are released. These electrons pass through the external circuit 
and reach the cathode electrolyte layer, and thus the current 
flows though the closed circuit. The reactions inside the FC 
stack can be summarized as follows: 


Reforming: CH4+H20 < CO + 3H2 (2a) 
Shifting: CO+H20 <> CO2 +H? (2b) 
Electrochemical: H2+1/2 O2 > H20 (2c) 


In the anodic recirculation SOFC system, STCR, defined as 


(3) 


NH0 


STCR = ——>——_ 
nco + NCH, 


is a very important parameter to evaluate carbon deposition in 
the reformer and FC stack. In terms of ejectors, the STCR can 
be evaluated by the entrainment ratio w, defined as [10] 


ms 
o = — (4) 

mp 
The entrainment ratio mainly depends on the ejector’s geome- 
tries and the operation conditions such as the primary flow 
pressure, secondary flow pressure and the back pressure 
(reformer pressure). For a given ejector, the behavior of entrain- 
ment ratio is strongly influenced by the back pressure. As shown 
in Fig. 2 [11,15], there are three different operational modes for 
different back pressure: the critical mode, the sub-critical mode 
and the back flow mode. The maximum entrainment ratio is 
reached in the critical mode and is independent of the back pres- 
sure. As the back pressure increases from its critical value, the 
entrainment ratio drops rapidly, and then becomes zero or nega- 
tive which implies that the ejector works at the back flow mode. 
From Fig. 2, it can also be seen that the entrainment ratio rises 


Subcritical 


ritical 
Critical Back flow 


mode | mode 


Primary pressure: 


P,<P2<P; 


Entrainment ratio 


. Back pressure 
Critical back pressure 


Fig. 2. Operational modes of ejector. 


as the primary flow pressure decreases. Since an ejector should 
normally work at the critical mode, ejector performance only at 
the critical mode operation will be discussed in this work. 


3. Governing equations 


A schematic view of the flow characteristic and pressure 
distribution along the flow direction in a typical critical mode 
working ejector are shown in Fig. 3 . The value of pressure is 
indicated by the degree of colori.e. deeper color represents larger 
pressure. The primary flow shocks at the converging—diverging 
nozzle throat (Section 1) and discharges from the nozzle exit 
supersonically, while its pressure is substantially dropped. The 
primary flow then expands after the nozzle exit introduces a 
series of oblique shocks in the suction chamber and accelerates 
the entrained secondary flow to choking condition at Section 3. 
The secondary flow mixes with the primary flow after shocking 
in the mixing chamber. The mixed flow will shock at the end of 
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Fig. 3. Schematic diagram showing the flow characteristic and pressure distri- 
bution in ejector. 
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the mixing chamber and in the diffuser to convert kinetic energy 
into pressure to reach a higher value than that of the secondary 
flow inlet one [16,17]. 

Without loss of generality, the following assumptions are 
made in developing an ejector model for SOFC operation: 


1. The velocity of the secondary flow inside the ejector is non- 
uniformly distributed in the radial direction and there exists 
velocity boundary layer near the inner walls of the ejector. 

2. The primary flow velocity is uniform in the radial direction, 
as there is no boundary layer and that the flow area of the 
primary flow is much smaller than that of the secondary flow. 

3. The primary flow is fully heated to the temperature of the 
secondary flow and the lost heat energy of secondary flow is 
negligible (i.e. Ts 3 = Ts,0; Tp3 = Ts,o), since the amount of 
the secondary flow is generally about five times that of the 
primary flow and the temperature of the secondary flow is 
much higher than that of the primary flow. 

4. Both the primary and the secondary flows are ideal gas inside 
adiabatic ejector walls. 

5. Pressure and temperature of both the primary and the sec- 
ondary flows are uniformly distributed inside ejector. 

6. The isentropic relations hold in calculating friction losses. 


3.1. Primary flow 


Using the relations between velocity and Mach number at 
each section 


Vp, i = Mp, 1 y kRgTp,i, 


and the fundamental equation for ideal gas 


i= 1,2,3 (5) 


eco (6) 
P one 
the governing equations for the primary flow in the 


converging—diverging nozzle and the suction chamber are given 
as follows: 


e Relations between Section 0 and Section 1 
Using the isentropic flow laws and energy balance, and 
taking the flow friction loss into consideration by multiplying 
an isentropic coefficient Wp, we obtain the formula for the 
primary mass flow rate 


2 
mp,1 = y WekRg (= peoAr/ Tro (7) 


where the average gas constant and density of the fuel are 
defined as 


) (k+1)/Q(k—-1)) 


R ni 
Re = Rubins o (8a) 
Jins oMo' 
and 
P P nl “Mo! 
fen P,0 P0 ipo (8b) 


~ RgTe.o  Tpo Rud nbo 


for a mixture inlet fuel, respectively. The detailed derivation 
is given in Appendix A. 
e Relations between Section | and Section 3 
Assuming that the primary flow fully expands after the noz- 
zle exit into the suction chamber, and using the isentropic flow 
and energy balance laws for the primary flow from Section 1 
to Section 3, we can obtain: 


(9) 


2(Pp,o/ Ps,o)*—P/* — 2 
mS = 


and 


Vp,3 = Mp 3y kRgTs,0 (10) 


e Diameter of primary flow at Section 3 


To determine the flow area of the primary flow in Section 3, 
the mass and energy balance for the primary flow from Section 
1 to Section 3 is required. The actual expansion diameter of 
primary flow at Section 3, Dp3, is expressed by 


D | 1 (2+(k-1)Mĝ, 
Dea = k+l 


(k+1)/(4K-1) 
(1) 


where &exp is a coefficient accounting for the frictional loss due 
to the mixing of two flows. The details of deriving Eqs. (9) and 
(11) are presented in Appendix B. 


3.2. Secondary flow 


According to geometry features of a fuel ejector, flow area 
of the secondary flow is much greater than that of the primary 
flow. Due to the turbulence and viscosity of fluid, the velocities 
inside the ejector are non-uniformly distributed and there exists 
a velocity boundary layer in the secondary flow near the inner 
walls of the ejector as shown in Fig. 4 . In the traditional 1D 
model, both velocities of the primary flow and secondary flow 
are treated as constant in the radial direction. It is clear that this 
modeling method will cause large errors in modeling of fuel 
ejectors. 


Actual velocity distribution 


anaa Traditional 1-D approach 
saeg Present model 


i 
i 
i 
1 
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1 

i 


' Mixing layer (M=1) `^, 


Vp 


Fig. 4. Velocity distribution at the Section 3. 
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We now find out a simple yet accurate function to approxi- 
mate the actual velocities inside the ejector. Considering that the 
function should satisfy (1) boundary conditions: r=0; vr = Vp.3 
and r= R3; vr = 0; (2) velocity characteristics of turbulent flow 
inside a pipe, we propose the following governing equation for 
the velocities of primary flow and secondary flow in Section 3 


Vp.3 O<r< Rp3 


Uy = R3 — r\ ™ 12 
i wa( 7 ") Rp3 <r < R3 e 
3 


where ny is an exponent to be determined as follows. 

At critical mode, it can be reasonably assumed that only 
the layer between the primary flow and secondary flow is in 
the choking condition (M= 1) and the layer is very thin [12]. 
Therefore, from assumption 3, the velocity and the radius of the 
mixing layer are vr = \/kRgTs,9 and r=Rp3 = Dp3/2, respec- 
tively. Substituting these values and Eq. (10) into Eq. (12), we 
obtain 

ln(1 — Rp,3/R3) 
is oa —InMp 3 Cal 

The graph of v; versus r is given as the dashed curve in Fig. 4. 
Clearly, it is quite close to the actual velocity distribution. 

Based on this 2D expressional velocity of the secondary flow, 
we define the mean mass flow rate of secondary flow at Section 
3 as 


R3 
jee 1 Dv; dA (14) 
Rp3 


where p stands for the average density of secondary flow. Sub- 
stituting Eq. (12) into Eq. (14) leads to the following integral 
equation: 


R3 r 1/ny 
ms 3 = 2nVoap | (1 = ) rdr (15) 
Rp3 R3 


By evaluating the integral of Eq. (15), the mass flow rate of the 
secondary flow can finally be expressed as 


n, R3? w Rp3 (ny+1)/ny 
ny +1 R3 


nR (i Res) OTD ao 
2ny + 1 R3 


ms,3 = 27 Vp 3p 


where the average density of the secondary flow is given as 


Pso _ Pso E;n oMo' 
RgTs,0 Ts,0 Ru) ;ns o 


pa (17) 


3.3. Energy balance 


For an ideal gas, the energy balance of the primary and sec- 
ondary flow in the ejector can be described by 


S “mp oC} Tp,o + Som oC} Tso = X Omi CT, + Eoss 
i i i 
(18) 


where the energy loss, Ejoss, of primary flow and secondary flow 
in the ejector can be approximated as: 


Ejoss = (1 — Eexp)mp,oVp 3 (19) 


Remark. With this modeling method, the velocity of the major 
portion of the gas inside the fuel ejector, secondary flow, is accu- 
rately modeled by a 2D function. Because of this improvement, 
the modeling of complicated mixing process in the mixing cham- 
ber and diffusing in the diffuser is avoided in the present model. 
Note that the characteristics of the mixing process in the ejector 
have not been well understood yet. Therefore, the present model 
will result in less uncertainties or errors in the modeling of fuel 
ejectors in the SOFC system. Of course, this will simplify com- 
putation significantly. For example, for fuel ejector geometries 
design the model only consists of seven algebraic equations: 
Eqs. (3), (7), (9)—(1 1), (13) and (16). 


4. Design and simulation procedures 
4.1. Design procedure 


The whole ejector design procedure is given as follows: 


1. Specify the operating conditions, including: 
e fuel inlet mass flow rate mp; 
e pressure and temperature of fuel inlet Ppo and Tpo; 
e pressure and temperature of anodic recycle gas Ps 0; Ts,0; 
e chemical composition of fuel and anodic recycle gas; 
e STCR. 

2. Determine the mass flow rate of anodic recycle gas ms from 
Eq. (3). 

3. Determine the nozzle throat area A; from Eq. (7). 

4. Compute the velocity of primary flow at mixing chamber 
inlet Vp3 from Eq. (10) 

5. Compute the diameter of active flow area of primary flow 
at mixing chamber inlet Dp3 from Eq. (11). 

6. Determine the radius of mixing chamber R3 from Eqs. (13) 
and (16). 

The remaining steps are to follow the standard procedures 
in the design manuals [11,18], which are listed below for 
the completeness of the paper. 

7. The angle employed in the divergent nozzle, a1, ranges from 
4° to 7.5°, with 5° to 6° most common. 

8. The distance of nozzle exit to mixing chamber inlet, Lg, is 
around 1.5D3 for best performance. The angle of suction 
chamber, a2, ranges from 3.5° to 5°. 

9. The mixing chamber is typically 3D3 to 5D3 long to accom- 
modate the shock pattern. 

10. The diffuser is always conical in shape with an included 
angle range of 5—12° and 8—10° is most common. Its length, 
La, is usually 4D3 to 12D3 in practice. 


A detail flowchart of the design procedure is shown in 
Fig. 5. 
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mp, STCR, Composition 


of fuel and anodic gas 
———>_ Equation (3) 
1 ms 


Ppo, Tpo, mp A 
Equation (7) 

Pso N 
——>| Equation (9) 


Tso 


——>| Equations (10), (11) 


Initial Value of D3 


—— | Equation (13) 


Equation (16) 


YES 4 D3 


Fig. 5. Calculation flowchart for fuel ejector geometries design. 


4.2. Performance simulation procedure 


For a given ejector geometry, its performance depends on 
the pressure Ppo and Ps,o, temperature Tpo and Ts.9, chemical 
composition of inlet fuel and anodic recycle gas, respectively. 
Since the main functions of the ejector in SOFC system are 
to meet the STCR at appropriate level and provide sufficient 
heat energy for the reforming reaction, the outputs at ejector 
exit are the entrainment ratio, STCR, temperature and chemical 
compositions. The procedure for performance simulation of the 
fuel ejector is very simple without iterative computation, which 
is given in the flowchart shown in Fig. 6. 


5. Model comparisons and validation 


The proposed models for ejector design and performance 
simulation will be validated from three aspects: design model 
validation, performance validation at on-design condition and 
performance validation at off-design conditions. 


5.1. Design model validation 


Table 1 shows typical values for an SOFC system [13]. Pure 
methane is used as fuel; the chemical compositions of the anodic 
exhaust are CH4, H20, CO2, CO and H2. Two coefficients, Wp 
and &exp, used in the design are both taken as 0.98. 

Starting with these data, the main parameters are obtained 
by following the proposed ejector design procedures as listed in 
Table 2 , in which the design data obtained by Marsano et al. 


Peo, Tpo, D : 
PERO ty Equation (7) 
Mp 
Pso } 
—| Equation (9) 
Mp3 
Tso : 
Equations (10), (11) 
Vp3, Dp3 
Ds y Equation (13) 


Ny 


Equation (16) 


Ms 3 


Equations (3) 


STCR 


Equations (18), (19) 


T4 


Fig. 6. Calculation flowchart for fuel ejector simulation. 


Composition of 


fuel and anodic gas 


[13] are also listed for comparison. It is shown from the table 
that the two critical ejector parameters: nozzle throat diameter 
and mixing chamber diameter agree fairly well. Even though the 
values of other parameters Da, Ls, Lm and Lg are different, they 
are not critical since 


1. given the same nozzle throat diameter and mixing chamber 
diameter, different design references gives different design 
parameters; 

2. different values of D4, Ls, Lm and Lg will not seriously affect 
on the ejector performance. 


Table 1 
Design point values of ejector 
Parameter Value 
Fuel inlet 
Composition (molar, %) CH; (100) 
Flow rate (kg s7!) 0.0094 
Pressure (bar) 10.06 
Temperature (K) 673 
Anodic recirculation gas 
Composition (molar, %) Hp (4.895) 
CO (3.785) 
H20 (61.74) 
CO2 (29.58) 
Pressure (bar) 3.8 
Temperature (K) 1280 
Steam to carbon ratio 2.4 
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Table 2 
Ejector design results 


Parameter Present model Marsano et al. [13] 


By calculation 


Nozzle, D; (mm) 3.31 3.31 

Mixing chamber, D3 (mm) 19.98 20.04 
By reference book 

D4 (mm) 53.5 100.8 

L; (mm) 40.0 — 

Lm (mm) 100.0 219.0 

La (mm) 239.8 459.5 


5.2. Performance validation at on-design condition 


For comparison purpose, ejector geometries and inputs for 
both the proposed and the Marsano et al.’s model are taken from 
Marsano et al.’s ejector [13]. The main results obtained for the 
ejector on-design operation conditions are listed in Table 3. It is 
observed that the primary mass flow rate, secondary mass flow 
rate, entrainment ratio and STCR for both models agree very 
well, where the maximum difference is only 1.94%. 

In the above two cases, the main advantage of the proposed 
model over Marsano’s is the reduction of computational cost. 
As stated in Section 3, only modeling the flow inside the suction 
chamber of ejector is required by the present model, while mod- 
eling of the other two chambers such as mixing chamber and 
diffuser is not needed. In contrast to the present model, tradi- 
tional 1D models need to apply the mass, energy and momentum 
conservation equations into each chamber to calculate the ejector 
performance. 


5.3. Performance validation at off-design condition 


Since the SOFC system load is usually adjusted through the 
fuel inlet flow rate, performance validation at off-design condi- 
tion can be established by varying the fuel inlet pressure. The 
general performance of a refrigeration ejector obtained by exper- 
iment [15] and computational fluid dynamics (CFD) simulations 
[16] is shown in Fig. 7 , where the primary mass flow rate mp 
rises linearly with the primary flow pressure in all the three oper- 
ating modes. In contrast, the behavior of secondary flow is more 
complex. 


e At low primary flow pressure, no secondary flow is entrained 
into the ejector. 


Table 3 
Ejector simulation results 


| | 
Back ISubcritical 
flow | mode | 


Critical mode 
(Design operation range) et 


- 
- 


Ejector primary pressure 


Critical primary pressure : 
(Fuel inlet pressure) 


Fig. 7. Ejector performance at different primary pressures. 


e As the primary flow pressure increases to the sub-critical 
region, the secondary flow is entrained proportionally and 
very sensible to the primary flow pressure. A maximum of 
secondary flow rate occurs when working conditions reaches 
the critical mode. 

e In the critical mode, the secondary flow rate decreases and 
then remains quite constantly in the high primary flow pres- 
sure regions. 


Since an ejector should normally work at the critical mode, 
ejector performance only at the critical mode operation will be 
discussed in the squeal. 

Due to lack of experimental data, the validation will only 
be carried out through comparing the results of the proposed 
model with the experimental data obtained with Freon as the 
working fluid, where the mass flow rate and entrainment ratio 
versus the primary flow pressure are shown in Fig. 7. Although 
the working fluid of the experimental data is not methane, it can 
still be treated as ideal gas such that the basic assumptions for the 
model are still valid. Therefore, the trends of the mass flow rate 
and entrainment ratio versus the primary flow pressure obtained 
from the fuel ejector are the same as that in Fig. 7. 

Keeping the fuel inlet (primary flow) temperature and anodic 
exhaust (secondary flow) conditions constant, the detail relation- 
ships between the two mass flow rates, entrainment ratio w and 
the fuel inlet pressure are shown in Fig. 8 . The mass flow rate 
of the primary flow increases linearly with the fuel inlet pres- 
sure in accordance with Eq. (7), while the secondary mass flow 
rate and also the slope of the curve decrease gradually. In the 


Marsano et al. [13] A =(Present’s — Marsano’s)/Marsano’s 


Parameter Present model 
Primary mass flow rate, mp (kg s~!) 0.0094 
Secondary mass flow rate, ms (kg s~!) 0.0689 


Entrainment ratio, w 7.34 
STCR 2.44 


0.0094 0 

0.0677 0.0177 
7.20 0.0194 
2.40 0.0167 
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Fig. 8. Entrainment ratio and mass flow rate at different fuel inlet pressures. 


high primary flow pressure region, the secondary mass flow rate 
is almost constant. These clearly agree with the critical mode 
experimental data shown in Fig. 7 well. 

However, in Fig. 6 of ref. [13], both the secondary mass 
flow rate and entrainment ratio increase with increasing the pri- 
mary flow pressure. Compared this to the experimental results, 
it is shown that the calculation error in Marsano’s ejector model 
enlarges as the primary flow pressure increases. Inaccuracy of 
the 1D model is much larger when applying it in off-design 
ejector performance evaluation than in on-design performance 
evaluation, mainly because the velocity distribution varies a lot 
in the off-design conditions. The 1D model suffers from the 
assumption of the uniform velocity distribution which makes it 
generate large error in the fuel ejector performance evaluation 
at off-design conditions. 

To explain in depth the phenomenon taking place inside the 
ejector, the relations of flow diameter Dp3 and velocity Vp3 of 
the primary flow at the cross Section 3 against the primary flow 
pressure are mapped in Fig. 9 . It shows that both the primary flow 
area and velocity increase as primary flow pressure increases. 
Consequently, the primary mass flow rate increases (Fig. 8). By 
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Fig. 9. Primary fluid velocity and flow diameter at Section 3 in different fuel 
inlet pressures. 


noting that the diameter of mixing chamber is fixed, the flow 
area for the secondary flow is decreased while Dp3 increases. 
Since the velocity of secondary flow is almost constant at critical 
working conditions, the secondary mass flow rate will decrease 
when primary flow pressure increases (Fig. 8). 


6. Performance influence of other parameters 


The ejector’s performance is also affected by temperature of 
primary flow and secondary flow, pressure of secondary flow and 
chemical composition. In this section, we will study the effects 
of these parameters on the ejector performance on flow rate of 
anodic recycle gas, STCR, pressure, temperature and chemical 
composition at ejector outlet. 


6.1. Influence of temperature 


Fig. 10 shows that the behavior of two mass flow rates 
and entrainment ratio at different fuel inlet temperature. It is 
observed that the primary mass flow rate decreases for increased 
primary temperature, while the primary temperature does not 
influence on the secondary mass flow rate. This is in line with 
ESDU [19] that the temperature of super heated primary flow 
does not affect the secondary mass flow rate. 

The relations between mass flow rates and secondary flow 
temperature are shown in Fig. 11 . It is seen that the sec- 
ondary mass flow rate decreases as its temperature increases. 
This can be explained as follows: (1) higher temperature 
means lower density for a fixed pressure gas, resulting in 
less entrained secondary mass flow; (2) gas dynamic viscosity 
increases when temperature increases and pressure is constant. 
According to fluid dynamics principles, it is known that mass 
flow rate is reciprocal to the dynamic viscosity in turbulent 
flow. 


6.2. Influence of secondary flow pressure 


The relations between mass flow rate and secondary flow 
pressure are shown in Fig. 12 , where the secondary flow 
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Fig. 10. Entrainment ratio and mass flow rate at different fuel inlet temperatures. 
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Fig. 11. Entrainment ratio and mass flow rate at different secondary tempera- 
tures. 


pressure has a strong influence on the secondary mass flow 
rate, as its density increases along with pressure when tem- 
perature is constant. Note that the fuel ejector is always 
working at its critical range (Fig. 7), i.e. the ejector exit pres- 
sure is always 1.015Ps when varying the secondary flow 
pressure. 

Fig. 13 shows the influence of both primary flow pressure 
and secondary flow pressure on the ejector performance. The 
STCR decreases as the primary flow pressure increases, but 
increases with the increasing secondary flow pressure. For a 
designed value of the STCR between 1.8 and 2.4, it can easily 
drop below the limit if primary flow pressure is high, with an 
evident risk of carbon deposition. 

Fig. 14 shows the influence of primary and secondary flow 
pressure on the ejector outlet temperature. The exit temperature 
decreases when primary flow pressure increases and secondary 
flow pressure decreases, as the temperature of primary flow is 
low while the secondary flow is high. 
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Fig. 12. Entrainment ratio and mass flow rate at different secondary pressures. 
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Fig. 13. STCR vs. fuel inlet pressure at different secondary pressures. 


6.3. Influence of secondary flow chemical composition 


The chemical composition of secondary flow depends on the 
fuel cell operating conditions. To analyze the effects of sec- 
ondary flow composition on ejector performances and STCR, 
we define the molar rate of H2O and CO7 at anodic recycle gas, 
Mpc, as: 


NH0 
NCO, 


Myc = (20) 


while the other two chemical compositions, Hz and CO, are con- 
stants of 4.895 and 3.785%, respectively. From Fig. 15 , it can be 
seen that the entrainment ratio w decreases with increased Myc, 
since the average density decreases as Myc increases because the 
molecular weight of H20 is lower than CO2. The ratio STCR 
increases with the increasing Myc but the slope of the curve 
decreases slightly and near constant at high Myc, in line with 
the trends of entrainment ratio and secondary flow composition. 
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Fig. 14. Ejector exit temperature vs. fuel inlet pressure at different secondary 
pressures. 
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Fig. 15. w and STCR in different chemical composition of secondary flow. 


7. Effects on SOFC system performance 


To study the influence of the fuel ejector on the anodic recir- 
culation SOFC system, typical SOFC system design values are 
listed in Table 4 . Starting with these data, the main task is to 
predict the performance of the SOFC system (entrainment ratio, 
STCR, temperature at each components, cell voltage and cell 
power) on varying fuel inlet pressure, cell pressure and fuel 
utilization. 


7.1. SOFC model 


The SOFC model developed by [20] is adopted in this paper. 
The model allows the evaluation of the design and off-design 
behavior of the SOFC system under the influences of the cell 
operation pressure and temperature, the fuel and air utilizations, 
the current density, voltage, etc. During the calculations, the 
main assumptions are listed as follows: 


1. Temperature within all the components of SOFC system is 
uniformly distributed. 

2. Cathode flow is composed of 21% O and 79% No. Fuel is 
pure CH4. 


Table 4 

Design point values of SOFC 

Parameter Value 

SOFC parameters 
Cathode thickness (cm) 0.035 
Electrolyte thickness (cm) 0.017 
Anode thickness (cm) 0.030 
Overall cell area (m°) 95 


Operation conditions 


Fuel utilization 0.85 
FC pressure (bar) 3.80 
FC pressure loss (kPa) Dal: 

Air inlet pressure (bar) 3.84 
Air inlet temperature (K) 1000 
Air flow rate (kg s7!) 0.47 


3. The reforming and shifting reactions are at equilibrium in the 
reformer and FC stack. 

4. Temperature of the gases at the outlet of the reformer and 
FC stack are equal to the reformer and FC stack temperature, 
respectively. 

5. The pressure loss in FC stack is equal to 1.5% of the FC 
operation pressure. 

6. The concentration loss is fixed at 1.5 x 1077 Q. 


In the simulation, we define the fuel utilization coefficient Us as 
sumed 
ima 
Ur = in in in (21) 
ny, + nco + ANH, 


consumed 


where ny, is the reaction rate of H2 in the FC stack and 


nity no and neu, are the molar flow rate of H2, CO and CH4 
into the FC stack, respectively. 

The simulation studies on the anodic recirculation SOFC sys- 
tem are carried out at different operation conditions by varying 
some key parameters, while keeping the cathode air pressure 
and temperature constant. The system consists of three model 
blocks: the fuel ejector, the reformer and the fuel cell stack. The 
parameters in the blocks are coupled and iterative computation 
is required: 


e The inputs of fuel ejector are fuel inlet pressure and temper- 
ature, anodic recycle gas pressure, temperature and chemical 
composition; the outputs are w, STCR, pressure, temperature 
and chemical composition at ejector exit. 

e The temperature and chemical composition of the FC inlet 
gas are determined by reactions taking place in the reformer. 

e In the FC stack, the values to be calculated are anode 
exhaust temperature, pressure, chemical composition, cath- 
ode exhaust temperature, chemical composition, FC voltage, 
current density, power, etc. 


7.2. Ejector performance: integrated versus stand-alone 


The ratios w and STCR by varying the fuel inlet pressure are 
shown in Fig. 16 . They exhibit different behaviors but have sim- 
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Fig. 16. w and STCR at different primary pressures. 
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Fig. 17. SOFC system temperature at different primary pressures. 


ilar trends with those obtained from Figs. 8 and 13, respectively. 
The differences are caused by that the temperature and chemi- 
cal composition of the anodic recycle gas are always fixed at the 
design points in the stand-alone ejector, while they are affected 
by the FC stack when the ejector is used as a component in SOFC 
system. Compared with the ratio w in Figs. 16 and 8, the max- 
imum deviation is only 0.72% at Ppo = 8.06 bar, which implies 
that the temperature and chemical composition of anodic recy- 
cle gas do not seriously affect the ejector entrainment ratio w. 
However they have significant effect on STCR; the maximum 
deviation between integrated and stand-alone ejectors is 18.8%. 


7.3. Components temperature 

The gas temperature at each component for different fuel 
inlet pressure is shown in Fig. 17 . In all cases, the tempera- 
ture decreases with increasing fuel inlet pressure, because of the 
primary flow is a low temperature fluid. 


7.4. FC performance 


The voltage, power and current density of fuel cell stack are 
obtained by varying fuel inlet pressure as shown in Fig. 18 . It is 
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Fig. 18. SOFC voltage, current density and power at different primary pressures. 
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Fig. 19. STCR vs. primary pressure at different secondary pressures. 


observed that the current density decreases with increasing fuel 
inlet pressure, which can be analyzed through the ejector off- 
design performance characteristic. On one hand, the flow rate of 
CH; into FC, ne,» increases as the fuel inlet pressure increases 
as indicated in Fig. 8; On the other hand, the secondary flow rate 
significant decreases as fuel inlet pressure increases from Fig. 8, 
which decreases the flow rate of nit, and ere It is found that 


the reduction rate in nin + nee, is higher than the increment rate 


in 4néy, so that ngensumed will finally decrease according to Eq. 


(21). This will cause the current density decreases as it is directly 
correlated with co aa From Fig. 18, it is also seen that the 
FC voltage rises with the fuel inlet pressure as decreasing the 
FC stack temperature leads to a lower cell open circuit potential. 
The FC power increases with increasing the fuel inlet pressure 
and is almost constant at high fuel pressure. 


7.5. STCR 


The STCR with the variation of the FC operation pressure is 
shown in Fig. 19 . It is quite different from that obtained in the 
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Fig. 20. STCR vs. primary pressure at different fuel utilizations. 
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stand-alone ejector (Fig. 13), especially in the ranges of low fuel 
inlet pressure and high FC pressure. This confirms that, for a FC 
pressure lower than 3.4 bar, the STCR can easily drop below the 
limited value and thus suffers from carbon deposition. 

The results obtained by varying the fuel utilization factor at 
different fuel inlet pressures are shown in Fig. 20. This indicates 
that the STCR is strongly influenced by Us. Note that Us also 
has a direct effect on conversion rate of Hz in FC stack, which 
affects the FC performances such as temperature, pressure and 
chemical composition at the anodic recycle gas. It shows that, 
for Uç lower than 0.7 the STCR can easily drop below the limited 
value. 


8. Conclusions 


In this paper, a new model for fuel ejector with large entrain- 
ment ratio, overheated working gas and low pressure increment 
in the anodic recirculation SOFC system is developed. The new 
model has no more than nine algebraic equations and can be 
applied in both ejector design and performance simulation. In 
contrast to the traditional 1D models, the new model has the 
following advantages: 


1. The new model uses a 2D function to approximate the real 
velocity distribution inside the ejector. This makes it much 
more accurate than the traditional 1D model in describing 
the flow field inside the ejector (Fig. 4). 

2. The modeling of flows in the mixing chamber and diffuser 
is not required by this model. In addition, throughout the 
modeling procedure, the energy conservation equation is not 
required for calculating the secondary mass flow rate. Both 
factors make present model simple and easy to apply. 

3. It can be used for not only the on-design case but also off- 
design performance evaluations for ejectors in the anodic 
recirculation SOFC system, while the traditional 1D model is 
incapable for the off-design ejector performance prediction. 


Furthermore, the ejector in the cases of on-design, off-design 
and integrated into an SOFC system performances have also 
been investigated. The main conclusions from the simulation 
results are summarized below: 


1. Mass flow rate of the fuel increases with increased the fuel 
inlet pressure. 

2. The behavior of STCR is similar to the secondary flow rate, 
which decreases as the fuel inlet pressure increases. 

3. The temperature and chemical composition of anodic recycle 
gas do not seriously affect the ejector entrainment ratio w, but 
have a significant impact on the ratio STCR. 

4. The current density of SOFC stack decreases but its power 
output increases with increasing the fuel inlet pressure at 
constant fuel utilization. 

5. The STCR can easily drop below the limited value causing 
the carbon deposition in high fuel inlet pressure regions. 


The proposed model is expected to have wide applications in 
geometries design and performance evaluation of the ejectors in 


SOFC systems. Utilizing the model, our immediate future work 
will focus on validating the model effectiveness on performance 
evaluation at off-design operation conditions by experimental 
studies and investigating the influence of the cathode side and 
the fuel conditions on the performance of the fuel ejector and 
SOFC systems. 


Appendix A. Primary flow relations between Section 0 
and Section 1 


Using the isentropic flow laws, the temperature and densities 
and pressures of primary flow at the Section 0 and Section 1, 
Tpo, Tp.1, Ppo, PPA» Pro, Pp, are related by 


Tp me a (A1) 
Tp,0 PPO 

and 

Tei (ay (A2) 
TP,0 Peo 

respectively. 


Since the velocity at Section 0 is negligible compared with 
the velocity at Section 1, energy conservation equation of fuel 
from Section 0 to Section 1, and the mass flow rate at Section 1 
are 


7 Ve 
CpTe.o = CpTp,1 + 2 (A3) 
and 
mp,1 = pp,1AtVp,1 (A4) 


respectively. 
Considering the Mach number is 1 at Section 1, invoking Eqs. 
(5), (6) and (A3), we obtain 


Teo k+1 
Eig ee Se A5 
Tp, 2 ve 
Substituting Eq. (A5) into Eq. (5), results 
OR Tay \ 1⁄7 
Vp = (—2 = A6 
P, 1 ( 4I ) (A6) 


Combining Eqs. (A1) and (A4)-(A6), and taking account the 
flow friction loss of an isentropic coefficient, Wp, we finally 
obtain Eq. (7) for the primary flow mass flow rate. 


Appendix B. Primary flow relations between Section 1 
and Section 3 


Since the primary flow expands fully in the suction chamber, 
ambient pressure of the expansion flow can be represented by 
the pressure of surrounding secondary flow Ps 9. Using the isen- 
tropic flow and energy balance laws for the primary flow from 
Section 1 to Section 3, we have: 


Te3 7 ew 


(B1) 
Tp,1 Ppt 
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Vou Vp.3 
CpTP,1 F D = CpTp3 + D (B2) 
where Tp3 and Vp3 are the temperature and velocity of primary 
flow at Section 3, respectively. 
Multiplying Eq. (B1) by (A2), and adding Eq. (B2) to (A3), 
we obtain 


T3 _ (Rey (B3) 
Tp,o Pp.o 

= VE 
CpTp.o = CpTp,3 + a (B4) 


Substituting Eqs. (5) and (6) into Eqs. (B3) and (B4) for Section 
3, and using assumption 3, we obtain Eqs. (9) and (10). 

To determine the flow area of the primary flow in the Section 
3, the mass balances for the primary flow from Section 1 to 
Section 3 is required which can be expressed: 


pp,1 At Vp.1 = pp,3Ap,3 Vp,3 (B5) 
using the isentropic flow laws, we obtain 

T: 1/(k-1) 
PRS _ (7) B6) 
pP,1 Tp,1 


Substituting Eqs. (5) and (6) into Eq. (B2) for Sections 1 and 3, 
results 


Vp3 _ Mp 3y kRgTp,3 B7 
Vp, VSkRegTP 1 
Tp _ 2+(k—1)Mp., 
T3  2+(k-!1) 


then substituted Eqs. (B6)—(B8) into Eq. (B5), the actual expan- 
sion diameter of primary flow at Section 3, Dp3, corrected by a 
coefficient is obtained and expressed in Eq. (11). 


(B8) 
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